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To date, a range of ion channels have been identified in 
chondrocytes using a number of different techniques, pre- 
dominantly electrophysiological and/or biomolecular; each 
of these has its advantages and disadvantages. Here we aim 
to compare and contrast the data available from biophysical 
and microarray experiments. This letter analyses recent tran- 
scriptomics datasets from chondrocytes, accessible from the 
European Bioinformatics Institute (EBI). We discuss whether 
such bioinformatic analysis of microarray datasets can poten- 
tially accelerate identification and discovery of ion channels in 
chondrocytes. The ion channels which appear most frequently 
across these microarray datasets are discussed, along with 
their possible functions. We discuss whether functional or pro- 
tein data exist which support the microarray data. A microarray 
experiment comparing gene expression in osteoarthritis and 
healthy cartilage is also analyzed and we verify the differen- 
tial expression of 2 of these genes experimentally, namely the 
genes encoding large calcium-activated potassium (BK) and 
aquaporin channels. 



Chondrocytes are the resident cells of cartilage which pro- 
duce, maintain, and degrade the extracellular matrix (ECM). 
Whilst the cells are non-excitable, they have already been shown 
to express a rich complement of ion channels. 1 Chondrocyte 
ion channels are involved in several critical functions including 
mechanotransduction 2 and apoptosis. 3 Many other functions 
in chondrocytes, and their precursors, have all been shown to 
directly involve ion channels. 4 

Expression data to date have largely been obtained using tra- 
ditional biochemical and physiological techniques such as flux 
studies, electrophysiology, and immunohistochemistry. More 
recently reverse transcription and real time (quantitative) poly- 
merase chain reaction (PCR) methods have also identified novel 
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ion channels in chondrocytes. 5 Over the past few years, microar- 
ray analysis or "transcriptomics" has added a new dimension to 
the study of gene expression. Whereas, in the past, one studied 
ion channels on a cell-by-cell, or antibody-by-antibody basis, 
this new "transcriptomic" or "expression profiling" technology 
potentially allows one to examine thousands of transcripts simul- 
taneously. The technique is clearly very powerful but there are 
drawbacks too, since every scientific technique has its limita- 
tions. However, by combining modern omics strategies and con- 
ventional biomolecular techniques one can gain a more thorough 
and subtle understanding of chondrocyte biology. 

The most widely used method for functional study of ion 
channels is patch-clamp electrophysiology. However, despite the 
prominence, historical dominance and power of this approach, 
it is a very slow method and in real cell systems (as opposed to 
heterologous expression systems) definitive identification of a 
specific ion channel is difficult. The identification of ion chan- 
nels by patch clamp electrophysiology usually depends on the 
availability of a specific pharmacological ligand or a very distinct 
biophysical profile. Both of these are actually far rarer than one 
may think. For example, some of the most important channels 
for the control of chondrocyte function are the transient receptor 
potential (TRP) channels. 6 These are difficult to identify elec- 
trophysiologically because they show weak voltage sensitivity and 
there are relatively few selective ligands. 

Traditionally, one has backed up the electrophysiology with 
immunohistochemistry or PCR studies, but even the strongest 
advocates of this approach would have to concede it is a slow 
process. Expression profiling with a microarray speeds up the 
characterization of ion channel transcription but, in addition to 
concerns over probe specificity, statistical validity and detection 
accuracy, there is a more fundamental issue. Even if particular 
ion channels are detected, and one accepts that this detection 
approach is valid, alone it gives no information on the role of that 
channel. Does the channel control membrane potential, volume, 
secretion, or have little role in that cell at all? The 2 data sets 
could be complementary to each other, however, if one makes 
the assumption that both experimental techniques are equally 
valid. The microarray data could direct electrophysiological 
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Table 1. Ion channels detected in all 10 of the microarray studies considered in this report. 



Gene 


Score/10 


Gene product description 


CLIC1 


10 


p64/intracellular chloride channel 1 


CLIC4 


10 


p64/intracellular chloride channel 4 


CLNS1A 


10 


Chloride channel nucleotide-sensitive, P-glycoprotein, pCln. 


KCNMA1 


10 


Large calcium-activated potassium channel (BK) 


SCN1B 


10 


Voltage-gated sodium channel, (3-subunit. Modulates activity of the voltage-gated sodium channel. 


VDAC1 


10 


voltage-dependent anion channel 1 


VDAC2 


10 


voltage-dependent anion channel 2 



experiments, by generating targets for study, and these biophysi- 
cal studies could examine the role of the expressed genes. It is well 
known that there is a statistically significant, but weak, correla- 
tion between total RNA abundance and protein expression. This 
has been particularly well studied in yeast 7 and E. coli where, 
from gene to gene and cell to cell, ratio of mRNA copy number 
to protein number can vary from 1:100 to 1:10000. 8 The degree 
of correlation also depends on the gene ontology and can be as 
high as R 2 =Q.5 or so. 9,10 In the following text, we will compare the 
results of recent transcriptomic studies with those of traditional 
methods and ask; do the data match each other and what does 
each dataset add to the other? 

In this letter we analyse and discuss the data from the fol- 
lowing 10 chondrocyte Affymetrix microarray datasets openly 
available on the European Bioinformatics Institute (EBI) expres- 
sion profiling database (EBI accession codes): E-GEOD-6119, 11 
E-GEOD-10024, 12 E-GEOD-10556, 13 E-GEOD-1277, 14 
E-GEOD-14402, 15 E-GEOD-16464, 16 E-GEOD-18052, 17 
E-GEOD-18394, 18 E-GEOD-7683, 19 E-GEOD-8077. 20 These 
datasets were found by searching for "chondrocyte (s)" OR "car- 
tilage and chondrocyte (s)". These datasets have been produced 
from expression profiling experiments of extracted or primary 
cultures of chondrocytes from "healthy" cartilage from a range 
of human, bovine, rat, and mouse tissue. 

A total of over 40,000 gene probe sets (reporters) have been 
used and these include those targeted at 200 or 300 ion chan- 
nel, or related, plasma membrane proteins. In this paper, we have 
accepted the author's detection limits and analysis where they 
have uploaded these to the European Bioinformatics Database. 
Where these are not given, we have processed the raw data with 
a Mas5 transform 21 using the open source Affymetrix Power 
Tools (Affymetrix). Probe annotations were derived from the 
Affymetrix annotations database files using custom Perl scripts 
(revisions "na31" except where stated). We have not analyzed 
data from other microarray platforms. 

Commonality between Datasets 

Typically, using standard assessment criteria, approximately 
50 to 100 channel targets are detected in chondrocytes. Not sur- 
prisingly the dataset from each study are rather different to each 
other, but there is a common set of 7 probes that are detected in 
all 10 of these studies. These are given in Table 1. It is important 
to note that this analysis does not select particular ion channel 



genes to study, but reveals all the genes commonly expressed in 
the above microarray datasets (note that there are over 300 ion 
channel gene probes encoded on each of the cited Affymetrix 
chips; accessory proteins such as the channel "tetramerization 
domain" proteins are excluded). With such a large number of ion 
channel transcripts on the Affymetrix chips one would expect the 
random co-detection of transcripts across datasets. To quantify 
this we used the following statistical analysis. The probability of 
more than n transcripts being co-localised in d datasets is given 
by: 

(Equation 1) 

where p is the probability of a given gene appearing in a data- 
sets, and there are a total of g genes on each array. This gives a 
p-value for 7 (or more) transcripts appearing in all 10 datasets by 
chance as p <le-l4. 

It should be noted that these microarray datasets were derived 
from different species (3 rat, 3 mouse, 3 human and 1 bovine) 
and there are potential differences in chondrocyte isolation pro- 
tocols. Constraining analysis to just rodent (6 datasets) returns 
a set of 23 commonly expressed ion channel genes (Table 2). 
Figure 1 quantitatively illustrates both the overlap of genes 
commonly expressed between species (Fig. 1A) and the over- 
lap between each of the transcripts from human, mouse and rat 
(Fig. IB, C, and D, respectively). It is evident that far more tran- 
scripts were detected in all 3 of the mouse datasets than in all 3 
of the human datasets. This could be for three reasons; firstly, it 
is possible that the sensitivity of the mouse chips is greater, but 
we have seen no specific evidence for this. Secondly, each of the 
protocols requires manual dissection and separation of chondro- 
cytes from the subchondral bone and adnexa. It is possible that 
mouse "chondrocyte" samples are inherently contaminated with 
non-chondrocyte tissue. In an electrophysiological or immuno- 
histochemical study such contamination would be relatively easy 
to detect, but in a biochemical protocol, where harvested tissue 
is macerated and then processed, it could be missed. Thirdly, 
it is conceivable that there are genuine phenotypic differences 
between chondrocytes in mice and other animals. Such differ- 
ences have been discussed elsewhere. 22,23 
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Figure 1. Commonality of gene transcript expression between datasets. (A) Commonality between species. These represent transcripts present in 
each of the 3 datasets for each of the human, mouse, and rat datasets (i.e., includes 9 datasets total). Far more are observed in rodent datasets (mouse 
especially) than human. (B) Commonality between human derived datasets. The human studies used tissue harvested from either adolescents receiv- 
ing limb length correction surgery (E-GEOD-1277), adults receiving ACL surgery (E-GEOD-16464), or post mortem (E-GEOD-10024). Samples were taken 
from knee (E-GEOD-16464), distal femur (E-GEOD-1277), or shoulder (E-GEOD-10024) and focussed on articular (E-GEOD-10024), mixed (E-GEOD-16464), 
or growth plate (E-GEOD-1277) chondrocytes. Chip Ids; E-GEOD-10024 used HG-U133A and E-GEOD-16464 used the slightly newer HG-U133A_Plus_2, 
but E-GEOD-1277 used the U95AV2 GeneChip. All 3 human studies used expanded chondrocytes, but E-GEOD-10024 and E-GEOD-16464 re-constituted 
those into 3D cultures. Extraction enzymes were collagenase P (E-GEOD-10024), clostridial collagenase and deoxyribonuclease I (E-GEOD-16464) and 
trypsin (E-GEOD-1277). (C) Commonality between mouse derived datasets. Rodent studies suffer from inherent difficulties in extraction of tissue since 
cartilage is thinner than larger animals. Tissue used from the microarray studies analysed in this letter came from a variety of joints from immature mice 
and are likely to include mixed chondrocyte phenotypes. Where stated explicitly, chondrocytes were expanded in monolayer cultures following col- 
lagenase based isolation (E-GEOD-8052 and E-GEOD-7683). All 3 studies (E-GEOD-10556, E-GEOD-18052 and E-GEOD-7683) used the same Affymetrix 
Mouse430_2 chips. (D) Commonality between rat derived datasets. The rat femoral head (E-GEOD-6119, E-GEOD-14402) or knee (E-GEOD-8077) tissue 
was harvested from a range of ages from one day old neonates (from which "only the outer two-thirds of cartilage" was used to select for articular chon- 
drocytes, E-GEOD-14402) to several month old rats (300-320g, E-GEOD-8077). Strain was either Wistar (E-GEOD-611 9), Sprague-Dawley (E-GEOD-8077) or 
not stated. E-GEOD-6119 and E-GEOD-14402 both used monolayer expanded chondrocytes following collagenase II based isolation. E-GEOD-6119 also 
included pronase, but E-GEOD-8077 used direct RNA extraction from macerated tissue. All the included rat studies used the Affymetrix Rat230_2 chips. 
One bovine dataset derived from chondrocytes 3D cultured from carpal bones of 3 to 6 mo old calves was also analyzed (E-GEOD-18394, Affymetrix 
Bovine chip, annotated with version "na29"). Since there was only 1 bovine chondrocyte dataset on EBI (albeit including a number of replicates) this is 
not included in the Venn diagrams. All other datasets were annotated with revision Affymetrix annotation version "na31". Note that each of the 3 species 
sets in (A) is equivalent to the commonly expressed regions of the Venn diagrams in (B, C, and D). 



The remainder of this letter focuses on whether the microar- 
ray data provides useful clues as to which channels are expressed 
in chondrocytes, and whether, therefore, initial processing of 
microarray data will improve the rate of channel detection in 
chondrocytes, or potentially other tissues. 

Chloride Channels 

The chloride channel superfamily' is huge and includes the 
large C1C family, CFTR protein, calcium-activated, volume- 
activated, P64 related chloride channels (Clns), and intracellular 



chloride channels (CLIC). In many systems chloride channels 
have been less well studied than cation channels, although 
they were some of the earliest ion channels identified in chon- 
drocytes. 24 ' 25 Analysis of the current datasets reveals that both 
p64-related (CLIO, CLIC4) and C1NS1A ion channels were 
detected in all 10 microarray experiments. The function of the 
corresponding channels are unknown, complex, controversial, 
or a combination of all these. CLIO, CLIC4, and other mem- 
bers of the CLIC family of proteins appear to be legitimate anion 
channels. 26 They are often referred to as "p64-related" simply 
because their earliest characterization appeared to be of a 64kDa 
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Figure 2. (see page 463) Immunohistochemical identification of 
KCNMA1 (BKct-subunit) and KCNB1 (BK pi-subunit) in sections of healthy 
and OA equine cartilage. The data from normal equines is reproduced 
with permission from Mobasheri et al. 45 Macroscopically normal articu- 
lar cartilage samples were obtained from weight-bearing regions of the 
metacarpophalangeal joints of horses of mixed breed, age, and sex. Joint 
tissues were sourced from an abattoir in Nantwich, Cheshire and Taunton 
Devon. Animals were euthanized for non-research purposes having 
been stunned before slaughter for meat in accordance with Welfare of 
Animals (Slaughter or Killing) Regulations 1995. Sections of normal (n=6) 
and OA (n=3) equine cartilage were probed for channel expression by 
immunohistochemistry essentially as previously described. 45 Sections 
were incubated overnight at 4°C with rabbit polyclonal antibodies to 
the KCNMA1 and KCNB1. Antibody dilutions used ranged from 1:200 
to 1:1500 in tris-buffered saline containing 1% bovine serum albumin. 
Slides were incubated with horseradish peroxidase-labelled polymer 
conjugated to affinity-purified goat anti-rabbit immunoglobulins. Cell 
nuclei were counterstained by incubation with aqueous haematoxylin 
(code no. S3309; Dako). Positive control samples were included from liver 
and kidney. Omission of primary antibody served as negative controls. 
Photomicrographs of immunostained tissue sections captured using 
Nikon Digital Sight DS-5M camera driven by Nikon Eclipsenet image 
capture software (Nikon). Positive staining is indicated by brown stain- 
ing and particular evident at middle/superficial zones. (C and D) Semi- 
quantified protein expression density. The largest increase in expression 
density (from data such as that illustrated in Fig. 2) is in the middle zone, 
for both KCNMA1 (C) and KCNMB1 (D) (a - and |3 - subunit respectively). 
Note that in OA tissue there was insufficient superficial data to quantify 
expression levels. 



Table 2. Ion channel gene IDs from all those gene transcripts detected in 
all 6 rodent microarray studies. 



Gene ID 


CACNA2D1 


KCNJ6 


CLCC1 


KCNK2 


CLCN3 


KCNK6 


CLCN4-2 


KCNMA1 


CLCN6 


SCN1B 


CLIC1 


TRPC1 


CLIC4 


TRPM7 


CLNS1A 


TRPV4 


KCNA6 


VDAC1 


KCNAB1 


VDAC2 


KCND1 


VDAC3 


KCNJ11 





Table 3. Ion channels detected in 9/10 of the microarray studies consid- 
ered in this report. 



Gene 


Score/1 0 


Gene product description 


VDAC3 


9 


voltage-dependent anion channel 3 


CLCN3, 7 


9 


chloride channel 3 and 7 


CLCC1 


9 


chloride channel 


TRPC1 


9 


transient receptor potential cation channel 


TRPV4 


9 


transient receptor potential cation channel 



The most striking observation is that a number of further chloride chan- 
nels were detected. Also, however, we find that 2 transient receptor poten- 
tial channels are commonly detected. TRPC1, a canonical TRP, and TRPV4, a 
vanilloid TRP channel. 



protein. 27 They appear to be of relatively low conductance for 
a chloride channel (8-40pS, depending on experimental con- 
ditions 28 " 30 ). Whilst, as their name implies, these channels can 
localize to intracellular compartments, 31,32 in some cell types 
they also appear to be in the plasma membrane and could serve 
a role in secretion. 32,33 One possibility is that the channel shuttles 
to and from the membrane in a cell cycle dependent way. 32 The 
pCln channel (sometimes referred to as P-glycoprotein, pICIn), 
also detected in 10 out of 10 datasets, was first identified by 
Paulmichl et al 34 as a putative rather ubiquitous volume-sensitive 
chloride channel. More recently, this volume regulator role has 
been refuted 35,36 and roles in gene-regulation and development 
have been proposed. 37,38 The controversy surrounding the nature 
of this channel is discussed in detail by Strange 39 and Furst. 40 The 
particular issue is that it cannot be clearly determined whether 
this is a volume-sensitive channel, another type of channel, or 
is a regulator of a channel endogenous in the various expression 
systems in which it has been studied. 39 Even the fundamental 
property of ion selectivity is controversial, since recent studies 
have shown a rather higher permeability of pCln to cations than 
would be expected for a chloride channel. 41 It is certainly pos- 
sible that one or other of these channels does constitute the chlo- 
ride channel identified by Tsuga et al 25 in chondrocytes, but this 



is unlikely since the channel observed by Tsuga et al. is more 
typical of a classical maxi-type chloride channel in electrophysi- 
ological and pharmacological terms. 26,42 

The Voltage Dependent Anion Channels 

The detection of both VDAC1 and VDAC2 is particularly 
interesting. The associated channels are thought of, generally, 
as mitochondrial ion channels, found in the outer membrane 
of this organelle. The proteins have also been detected in the 
plasma membrane, where they exhibit voltage-gated anion chan- 
nel activity (based on data in the NCBInr and UniProt data- 
bases). VDAC channels are also implicated in apoptosis and, as 
such, they will be of profound importance to all cells in which 
they are expressed. It has also now been suggested by a number 
of authors that some VDAC protein expression is also of plasma 
membrane ion channels. 43,44 The phenotype of this channel is 
that of a maxi-chloride channel. This would be very much in 
line with the original chondrocyte chloride channel work of 
Tsuga 25 and Sugimoto. 24 The channel identified in these stud- 
ies was the maxi-chloride channel, which is remarkably similar 
to the maxi-Cl/VDAC channel. In our own unpublished work, 
we see clear expression of a large conductance, niflumic acid 
and SITS-sensitive chloride channel which appears likely to be 
maxi-chloride. 

The Large Calcium-Activated Potassium Channel 

It is no surprise that KCNMA1 (BK) has been detected 
in all of the 10 datasets analyzed here. Currents either 
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Figure 2. For figure legend, see page 462. 
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Figure 3. Water Permeability (aquaporin) Assay. Water permeability can be calculated from the initial slope of the relative volume (V/VJ against time 
curve. Where V is the volume at time f and V is the volume at time zero. This is the accepted physiological assay for aquaporin expression. (A) Permeability 
is 30±3% (p<0.05, n=4) greater in chondrocytes from dogs with osteoarthritis (OA). TEA (a blocker of AQP1 5 " 4 ), bumetanide 54 {plC50 5.17+0.11uM, n=6, 
"Bumex", a blocker of AQP1 and 4), and mercuric chloride (HgCI 2 a non-specific AQP blocker reversed by 2-mercaptoethanol) ME) are included to deter- 
mine AQP type (B and C). Chondrocytes were harvested from canine clinical waste tissue with owner consent. Cells were placed in a "physiological 
saline" solution including 120 mM sucrose (osmolarity 300mOsm),then moved to an identical physiological saline without the sucrose. Cells at first swell 
as water enters the cell due to osmosis. Live cell imaging was achieved with a Nikon Diaphot microscope equipped with a Sony ICX098QB high sensitivity 
CCD. Images were analysed offline with Image! Volume was calculated from the 2D surface area (A) of the cell disc by assuming the cell is approximately 
spherical as described previously, 6 using the following equation: / p- -\ 3 

Vol = -TC 




(Eqn.2) 

Except where stated, data are presented normalized for starting volume (V0) as V/V0, where V is the volume at time t. Visual data were analyzed with 
ImageJ and ANOVA performed with SPSS (SPSS Inc.). Note that canine tissue was harvested from clinical waste tissue with Local Ethical Approval, no 
dogs were harmed for the study. 



Table 4. Ion channels and porins significantly changed in OA. 



Gene Symbol 


Encoded ion channel 


Abundance 
ratio 


p-value 


AQP1 


Aquaporin 1 


39.8 


7.3E-22 


KCNK5 


K2P5.1 (Task-2) 


-4.7 


4.8E-16 


KCNMA1 


KCa 1.1 (BK) 


3.1 


5.0E-10 


KCNN4 


KCa3.1 (IK) 


10.2 


2.0E-1 7 


KCNT2 


BK channel subunit (KCa4.2) 


-2.2 


2.0E-07 


SCNN1A 


ENaC 


-3.6 


1 .2E-08 


TMEM16A 


Calcium activated 
chloride channel 


3.2 


1 .4E-20 



Negative is fold decrease in abundance ratio, positive is fold increase in 
abundance ratio. Data from Karlsson et al 2010, 50 published in supplemen- 
tary materials. Note that whilst there were probes for TRPV4 present on the 
chip, the huge variability between abundance scores between samples for 
this particular probe set make it unlikely that any change would be detect- 
able, even if there was one. 



broadly, or specifically, identified as being through BK chan- 
nels have been described in a number of papers. In our own 
work we have shown not only that BK currents are present by 
electrophysiology, but also that the KCNMA1 was detectable 
by immunohistochemical methods. 45 Interestingly, whilst the 
BK function modifying (3-subunit KCNMB1 was detected 
in our own immunohistochemical studies it was not detected 
in any of the 10 microarray studies discussed here. There was 
detection of KCNMB2 and KCNB4 in 2 and 3 out of the 
10 studies respectively. Activation of BK results in such large 
currents that it is likely to be involved with regulation of intra- 
cellular osmolarity and volume. This is important, since it 
known to be activated by stretch. 45 Early work did not show 
whether this activation was direct or indirect via an increase of 
intracellular calcium ions, but BK is, so far, the only candidate 
for this role. 
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The Sodium Channel (3-Subunit 

The detection of the sodium (31-subunit SCN1B in all 10 
studies is interesting, since the a-subunits were undetectable in 
most studies. This subunit has been shown, in neurons to convey 
subtle changes to expression patterns and functional properties of 
voltage-gated sodium channels. 46 There has been 1 electrophysi- 
ological study, which identified voltage-gated sodium channels 
in chondrocytes. 24 This work has not been followed up on and, 
since chondrocytes do not fire action potentials, it is difficult to 
see what the role of a voltage-gated sodium channel might be. For 
completeness: SCN7A (Na 2.1), a somewhat atypical Na* chan- 
nel was detected in 4/10 experiments, and the classical sodium 
channel SCN2A1 in 3/10, and the transient type sodium channel 
SCN11A and SCN10A in 2/10 datasets. SCN5A and 3A were 
detected in 1/10. 

Transient Receptor Potential (TRP) Channels 

The next most commonly identified ion channel transcripts 
are given in Table 3. This list includes TRPC1, of the canoni- 
cal TRP, and TRPV4, a vanilloid channel, both of which were 
identified in 9/10 microarray studies. In traditional protein and 
pharmacological studies, however, 3 TRP channels have been 
identified in chondrocytes TRPV4, 4748 TRP V5, 6,48 and TRPV6 48 
but not TRPC1. We have not investigated TRPC1 ourselves, 
however. The TRP channels identified pharmacologically in 
chondrocytes have different, but related roles. TRPV5 appears 
important for setting the membrane potential, crucial to mainte- 
nance of cell volume in chondrocytes. 6 TRPV4, however, appears 
to be critical for allowing entry of Ca 2+ and activation of BK 
channels during imposed volume increase, and thus the process 
of regulatory volume decrease in epithelial cells. 49 

Changes of Transcript Levels with 
Onset of Osteoarthritis 

To use traditional (electrophysiological) methods to iden- 
tify changes in all channels expressed between populations of 
cells from different tissues is probably not feasible unless high 
throughput automated ion channel recording equipment can be 
adapted for this purpose. However, observation of changes in cel- 
lular properties and mathematical models may provide clues as 
to which channels have changed in OA. An alternative approach 
would be to use microarray comparisons between normal and 
OA tissue, detect changes in channel mRNA abundance and 
then follow-up with functional or histochemical experiments. 
Generally, differential expression patterns from microarray stud- 
ies are used to identify changed pathways, however, this process 
could miss changes which take place in individual ion channels 
not associated with an established pathway. The next microarray 
study considered in this article specifically compared transcript 
abundance from chondrocytes in human normal and OA carti- 
lage. 50 The authors focused on (and verified) changes in many 
classic cartilage phenotypic markers, but did not specifically con- 
sider changes in other proteins such as ion channels. Probing this 



list for ion channels and porins significantly changed (greater 
than 2-fold change) produces the data in Table 4. The data clearly 
show an approximately 3-fold change in ENaC, TMEM16A and 
BK (KCNMA1) transcript abundance and a 38-fold change in 
the aquaporin AQP1 channel transcript abundance. Interestingly 
each of these channels are important for chondrocyte cell volume 
control. 2 ' 55 From this, by way of proof of principle we verified 
BK channel changes with immunohistochemistry and aquaporin 
channel changes with a functional assay. 

We analysed, using immunohistochemistry, whether protein 
expression of BK is increased. Tissue was taken from stifle joints 
of horses with and without OA (see Fig. 2A and B for repre- 
sentative example). We investigated the expression of both a - 
and (3-subunits of BK (only KCNMA1 was included in Karlsson 
2010, 50 there were no probe sets for KCNMB1 on the chip). 
Semi-quantitative analysis of protein expression density shows 
that both BK subunits were significantly increased in OA, in the 
middle zone (Fig. 2C and D). 

The accepted assay for aquaporin expression is that developed 
by Preston et al." Cells are challenged with a hypotonic solu- 
tion, causing them to swell, and the rate of swell can be mea- 
sured to determine the aquaporin expression. The bioinformatics 
showed a 3 8 -fold increase in AQP1 transcript abundance in OA 
and we see a significant increase in functional aquaporin expres- 
sion in tissue from OA joints (Fig. 3A), although the increase is 
much smaller than the change in AQP1 transcript abundance. 
Current pharmacological tools do not allow categorical determi- 
nation of aquaporin subtype, however, here, water permeability 
of chondrocytes was blocked by concentrations of bumetanide, 
TEA, and HgCl, consistent with that expected for AQP1 53 ' 54 
(Fig. 3B and C). 

In summary, there appears to be considerable agreement 
between transcriptomic studies and physiological or immunohis- 
tochemical studies. It would seem that most channels common 
to all 10 datasets can be identified by these other techniques. 
There are examples, however, of proteins which have been iden- 
tified in chondrocytes yet show up in few datasets. For example, 
the ASIC channel (ACCN2 gene) has been shown by immuno- 
histochemistry and rt-PCR yet shows up in only one of the three 
rat datasets discussed here. 52 It is also notable that KATP potas- 
sium channels, previously detected both by immunohistochem- 
istry and electrophysiology were not detected in the microarray 
studies. 55,56 Therefore, combining these approaches should mas- 
sively speed up the rate of discovery of ion channels in cell types 
which can be isolated in sufficient quantities to perform such 
studies. There are tissues, such as the brain, where cell types are 
too intermingled to allow identification of unique cell types, but 
for many tissues in the musculoskeletal system (or cell lines), the 
combination of transcriptomics and protein studies seems ideal. 
With regard to OA, this strategy has allowed us to very rap- 
idly pinpoint some phenotypic changes in the expression of two 
important channels in OA: an aquaporin and the BK channel. 
Further protein and functional experiments will be necessary to 
establish whether the other K., channels are also altered, and in 
particular whether these changes contribute to or result from 
progression of OA. 
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